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Abstract. Some peculiarities of the exploitation of the entropy inequality in 
case of weakly nonlocal continuum theories are investigated and refined. As 
an example it is shown that the proper application of the Liu procedure leads 
to the Ginzburg-Landau equation in case of a weakly nonlocal extension of the 
constitutive space of the simplest internal variable theories. 



1. Introduction 

Ginzburg-Landau equation and its variants appear in different fields of physics 
and are applied to several phenomena. Its physical content and the way to obtain 
it is sound and transparent. The traditional derivation of the equation comprises 
two main ingredients (see e.g. pQ) 

- The static, equilibrium part of the equation is derived from a variational 
principle. 

- The dynamic part is added by stability arguments (relaxational form). 

The two parts are connected loosely and in an ad-hoc manner. As a classical 
field equation defined on nonrelativistic space-time, the Ginzburg-Landau equation 
should be compatible with the general balance and constitutive structure of contin- 
uum physics. Recently there has been several efforts to give a uniform reasoning of 
the equation on pure thermodynamic ground and to generalize the method of the 
derivation [21 El EH E|- The treatment of Ginzburg-Landau equation is a kind of test 
of weakly nonlocal (gradient) thermodynamic theories. 

In some previous works the possible role of internal, dynamic variables in weakly 
nonlocal extensions of the constitutive state space was investigated both with the 
heuristic method of classical irreversible thermodynamics and with the more exact 
Liu procedure 01 [5] . It was found that with the help of special dynamic variables 
(current multipliers) appearing in a generalized entropy current one can find an 
equation that is similar to the Ginzburg-Landau equation, but is not the same. 
The equation was called thermodynamic Ginzburg-Landau equation. The static 
part of the Ginzburg-Landau equation gives nonhomogeneous solutions in case of 
uniform boundary conditions, therefore it is a pattern forming equation. On the 
other hand, thermodynamic Ginzburg-Landau is not a pattern forming equation, 
its corresponding solutions are homogeneous. This is an important difference and 
corresponds well to the fact that with internal variables one introduces a local the- 
ory, all the functions in the basic state space are local. Therefore any nonlocalitics 
formulated by an internal variable theory are relocalized in this sense. Pattern 
forming theories, like those based on the Ginzburg-Landau equation, cannot be 
relocalized. 
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In this paper it is shown, that there is a non-equilibrium thermodynamic ap- 
proach also for pattern forming equations, however, it is not relocalizable, indepen- 
dent of any kind of internal variables. A pure thermodynamic derivation can unify 
the mentioned two parts of the derivation, there is no need to postulate any kind 
of variational principle. The Euler-Lagrange form of the static part of the equation 
turns out to be a consequence of the Second Law. 

In the paper the key ingredients of the methodology and the mathematical back- 
ground to exploit the Second Law inequality is shortly summarized in weakly non- 
local continuum theories (gradient theories, theories with coarse grained thermo- 
dynamic potentials, phase-field models, etc.). After the introduction, an overview 
of the methods to exploit the Second Law inequality is given. We call the attention 
to the essential differences between the local and weakly nonlocal theories and a 
constructive way to arrive solvable Liu equations is given. The most important dif- 
ference is that the weakly nonlocal extension of the constitutive state space, implies 
that some space derivatives of the constraints (evolution equations) are to be consid- 
ered as constraints, too. A simple example show the application of Liu procedure. 
In the following section Liu procedure is applied to derive the Ginzburg-Landau 
equation. Different kind of generalizations are treated, the difference between the 
relocalized and pattern forming (non relocalizabile) theories is formulated. A new, 
generalized formulation of Liu's theorem and its direct proof from Farkas' lemma 
is given in an Appendix. 

2. Methods to exploit the Second Law inequality - procedures of 
cit, coleman-noll and llu 

In every non-equilibrium thermodynamic theories an important theoretical prob- 
lem is to formulate the correct form of the evolution equations taking into account 
the requirement of the entropy inequality. The most predictive methodological solu- 
tion of the problem is far from being trivial and originated from Coleman and Mizel. 
They essentially reverse the way of thinking: one should look for the solution of the 
entropy inequality taking into account the evolution equations as constraints OH]- 
In continuum physics the dynamic equations are given in a determined form (e.g. as 
balances of extensives), except some constitutive, material functions. The task is to 
ensure the nonnegativity of the entropy production with appropriate constitutive 
assumptions. Therefore one should specify the undetermined material functions 
such, that in case of all possible solutions of the dynamic equations the form of 
the constitutive functions, the material properties, ensure the nonnegativity of the 
entropy production. In case of weakly nonlocal theories the entropy current plays 
a distinguished role. The entropy and the entropy current are both constitutive 
and are to be determined according to the above requirement (as it was suggested 
in extended rational thermodynamics 8 ). The entropy possibly should preserve 
its potential character in a general sense, therefore solving the above problem (e.g. 
Liu equations) the practical aim is a simplification in such a way that every consti- 
tutive quantity including the entropy current could be calculated from the entropy 
function. 

There are three basic methods to exploit the entropy inequality. 

- Heuristic. This is the force-current method of classical irreversible ther- 
modynamics (CIT). In classical problems, in case of simple state spaces, 
the form of the entropy production is quadratic and the inequality can be 
solved. The method can be justified by the Liu procedure in case of the 
traditional, simplest state spaces , and the method can be generalized to 
include non-classical entropy currents |5HTT)]. 
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- Coleman-Noll procedure. In the Coleman-Noll procedure one exploits the 
constraints (e.g. dynamic equations) directly, substituting them into the 
entropy inequality. The degenerate form of Liu's theorem ( H5.5|l is ap- 
plied. One usually assumes a specific form of the entropy current. There 
are essentially two choices here. The entropy current can be the classical 
(is = jq/r) , or a generalized one. In weakly nonlocal considerations both 
classical (see e.g. and generalized forms are applied. Generalizations 
of the entropy currents (or currents of other thermodynamic potentials) can 
be suggested on different grounds and they give good results with the pro- 
cedure huei . 

- Liu procedure. With Liu procedure one applies Liu's theorem with Lagrange- 
Farkas multipliers ( 1 15. 31 in the Appendix). At the first glance the applica- 
tion of this method seems to have only practical advantages. However, as 
the Lagrange multiplier method preserves the simple form of the constraints 
in question in extremum problems, with Lagrange- Farkas multipliers one 
can preserve and exploit the structure of the constraint and the entropy 
inequality. The question is not purely mathematical, because there are 
cases where the multipliers cannot be eliminated and they can get physi- 
cal significance. Moreover, an inevitable advantage of Liu's method is that 
the structure of entropy inequality makes possible to solve completely the 
physical problem. 

The train of thought is the following. The entropy current is considered 
as an independent constitutive quantity. With a proper choice of the con- 
stitutive space we can solve the Liu equations and determine the entropy 
current. Hence the entropy inequality simplifies considerably. The point of 
view of Onsagerian CIT is important here: with a proper identification of 
thermodynamic currents and forces the resulted entropy inequality can be 
solved, determining all constitutive quantities. This ensures, that our the- 
ory is independent of further artificial constraints, the entropy inequality 
becomes a consequence of material properties. 

In weakly nonlocal continuum physical calculations with Liu procedure one 
should consider some additional practical rules. There the constraints are (par- 
tial) differential equations. The functions in the differential equations form the 
basic state space. The constitutive quantities depend on these functions, on the 
basic state and on some of its derivatives. These derivatives are locally independent 
therefore the problem is algebraically manageable. The basic state variables and 
some of its derivatives can be included into the constitutive state space (or sim- 
ply state space |13|L into the domain of the constitutive functions. The entropy 
inequality with its special balance form determines the independent variables of 
the algebraic problem: those are the derivatives of the constitutive state, the so 
called process directions. The choice of the constitutive state space is crucial and 
determines the restricted constitutive functions, after applying Liu procedure, the 
peculiarity of weakly nonlocal theories is that depending on the particular state 
space, some space derivatives of the original constraints (e.g. dynamic equations) 
further restrict the process direction space, therefore they should be considered in 
Liu's theorem as additional constraints. In the following we will show some exam- 
ples to clarify the most important practical rules in the application of the formalism. 
One can find other examples on the application of derivative constraints in |14l 115) . 

Remark 2.1. The mentioned exploitation methods of the Second Law, being al- 
gebraic, are essentially independent of the solvability of the dynamic equations, 
whether the associated problems are well or ill posed, on the applied function spaces, 
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etc... For example the number of the equations can be less than the number of vari- 
ables in the basic state space. The wanted fields, the functions searched in the final 
resulting differential equation can be different from the basic state. 

Example 1. In this example the basic state space is formed by two times differen- 
tiate real functions x : M i— ► R. The constitutive space is spanned by the basic state 
and its derivative {x,x'). We are looking for scalar valued differ entiable functions 
F and S as lying in the constitutive space so that 

S'(x,x') > 

for all {x,x r ) satisfying the constraint 

(1) F{x,x') = 0. 

Evidently S'{x,x') = d\Sx' + d 2 Sx" , where d n denotes the partial derivative 
according to the n-th variable. Therefore, the space of the process directions (the 
space of independent variables in Liu's theorem) is spanned by x" . We are looking 
for conditions on S and F that the above inequality should be true for all (x,x r ) 
solving Q), but independently of the values of x" . The degenerate case of Liu's 
theorem gives some conditions. The single Liu equation is 

d 2 S = 0. 

Therefore S is independent on x' . The dissipation inequality can be written in 
the following simple form 

(2) dxSx' = ^-{x)x' > 

dx 

The above inequality does not give any condition for F. However, let us observe, 
that one of our previous assumptions was too strong. The process direction variable 
x" is not really independent on the state space, the derivative of gives a further 
restriction 

(3) diFx' + d 2 Fx" = 0. 

Considering this condition we apply Liu's theorem ( H5. S\) with the multiplier 
method, introducing the multipliers X\ and X 2 for the constraints £7J) and 0) re- 
spectively 

d x Sx' + d 2 Sx" - As^iiV + d 2 Fx") - XiF = 

= - A 2 5i^)x' + (d 2 S - X 2 d 2 F)x" - X X F > 

Therefore we can read the Liu equation as follows 

\ 2 d 2 F -0 2 S = 

Expressing the multiplier and substituting into the dissipation inequality we get 

d!Sx' - X 2 d!Fx' - X X F = (diS - d 2 S(d 2 F)- 1 d 1 F)x / — XiF > 

In this example we face to a partially degenerate case, hence with X± = we can 
give the general solution of the above inequality, as 

diS - d 2 S{d 2 F)- 1 d 1 F = L(x, x')x', 

where L is nonnegative. Given a function S we can calculate F, with appropriate 
conditions on L . For example if S{x,x') — x ■ x' and L — constant, thenF{x,x') — 
f(x L ~ 1 x') is a solution of the above equation for any / : I h> R. 
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3. Weakly nonlocal non-equilibrium thermodynamics - 
Ginzburg-Landau and thermodynamic Ginzburg-Landau equations 

Let us denote an internal variable (e.g. an order parameter of a second order 
phase transition) characterizing the microstructure of the material by £. In this 
case the Ginzburg-Landau equation can be written as 

(4) d t £ = -7ir e + 7 2 A£, 

where is the partial derivative of the appropriate thermodynamic potential (e.g. 

= where / is the free energy), and <9 t denotes the partial time derivative. 
Here we assumed that the internal variable £ is not related to the mechanical 
motion, therefore the choice of the frame (partial or substantial time derivatives) 
is irrelevant. 71 and 72 are material coefficients. The usual form of the Ginzburg- 
Landau free energy density is 

(5) /(£,V£) = /o(£)-7(V£)72, 

where 7 is a material coefficient, fo is the static (equilibrium) free energy and so 
L^ = /'(£). Gurtin gave a method to deduce equation from pure thermody- 
namic considerations, with the concept of microforce balance and showed that some 
additional terms should appear with a characteristic structure pi 116). The gener- 
alized form of the equation together with the characteristic term is the following 



(6) dtt = -71 r 5 + 72AC + 7aA3 t £, 

Such kind of terms appear in connection of several different phenomena and not 
only in case of the Ginzburg-Landau equation [171 118| . E. g. the Guyer-Krumhansl 
equation of heat conduction can be considered as a Cattaneo-Vernotte type wave 
heat conduction equation supplemented by a Gurtin term. 

It was argued that the Ginzburg-Landau equation is the first nonlocal extension 
of any kind of equation for an internal variable 01 H] and its characteristic func- 
tional form can be derived from the requirement of compatibility with the Second 
Law, without referring to variational principles. Hence, the reason of its wide- 
range applicability is well founded, because any internal variable that can charac- 
terize a material structure and is independent of other requirements should fulfill 
a Ginzburg-Landau equation in the first nonlocal approximation. The arguments 
were supported by calculations based on Liu's theorem. However, in an internal 
variable, completely relocalized theory one cannot derive directly Q), but only a 
very similar equation, that was called thermodynamic Ginzburg-Landau equation 

(7) <9 t £ = -7tiL ? + 7 T2 Arj, 

where 7ti and 7T2 are material coefficients. Moreover, the entropy (free energy) 
was proved to be gradient independent. One can see, that the equations {2J and 
(0 are similar but not the same at all. The essential qualitative difference is 
that static (equilibrium) solutions of the thermodynamic Ginzburg-Landau equa- 
tion with homogeneous boundary conditions are homogeneous but static solutions 
of the original Ginzburg-Landau equation with the same boundary conditions are 
not, they can form structures. The situation is well known and understood in super- 
conductors. The London equation (corresponding to the thermodynamic Ginzburg- 
Landau) does not determine the penetration length of the magnetic field, however, 
the Ginzburg-Landau equation gives that |19j . 

In the following, applying Liu's procedure with the methodology described in the 
previous section, we will derive the Ginzburg-Landau equations from very general 
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assumptions and show that the Gurtin terms are consequences of pure thermody- 
namic considerations, without referring new concepts, like the configurational force 
balance or virtual power, etc.. 

We are looking for a dynamic equation of £ in the following general form 

(8) d£-F= 0, 

where J 7 is a constitutive function, which form is to be restricted by the Second 
Law. The basic state space is spanned by £. Let us assume that the constitutive 
space is spanned by £, V£ and V 2 £. In this case the entropy inequality will be 

d t s + Vj s = 

dis dt$ + d 2 s ■ Vdt£ + d 3 s : V 2 d£ + dtf. • V£ + d 2 j s : V 2 £ + d 3 j s - : V 3 £ > 0. 

One can see, that the space of the process directions (independent variables) is 
spanned by <9t£, V(?t£, V 2 <9t£ and V 3 £. Moreover, let us observe that these variables 
are not really independent, the gradient of JHJ connect them. Therefore, in addition 
to JSJ one should consider the following constraint 

(9) Vd t £ + VJF = Vd t £ + dxTVS, + d 2 T ■ V 2 £ + d 3 T : V 3 £ = 0. 

Introducing T± and T 2 Lagrange-Farkas multipliers for the constraints JSJ and 
© respectively, one can get the following Liu equations 

dia = ri, 
d 2 s = r a , 

d 3 s = 0, 

(d 3 j s ~r 2 d 3 f) s = o. 

Here the superscript s denotes the symmetric part of the corresponding tensor. 
The first two equations determine the multipliers. From the third equation follows, 
that the entropy does not depend on the second derivative of £. Taking into account 
these requirements one can give a solution of the fourth equation and determine 
the entropy current as 

(10) v& v 2 o = d 2S (t, vo^(€, ve, v 2 o + j (f, V£). 

For the sake of clarity we explicitly denoted the variables of the corresponding 
functions. With the above solution of the Liu equations the dissipation inequality 
can be simplified considerably 

(11) V-j + (V-<9 2 s-dis)- J">0. 

Assuming, that jo = one can give the general solution of the above inequal- 
ity. That solution can be interpreted by the well known traditional method of 
irreversible thermodynamics, choosing appropriate forces and currents. Therefore, 
T is the constitutive quantity to be determined (thermodynamic current) and it 
should be proportional to the given one (force) 

(12) d t i = T = L(V ■ d 2 s - d lS ) 

with a nonnegative state dependent constitutive function L. s(£, V£) is a given 
entropy function (determined from static measurements). (11 2t is the Ginzburg- 
Landau equation, and one can get back the very traditional J3J form using the 
specific entropy functional with a form like JSJl and dealing with a strictly linear 
theory in a thermodynamic sense, where L is a constant function. The choice of 
the right thermodynamic potential (e.g. entropy or free energy) depends on the 
boundary conditions [Q. 
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If we do not restrict the space of the independent variables by JSJ, by the de- 
rivative of the original constraint, then after an easy calculations one can get the 
following form of the dissipation inequality 

(13) V • j (C, V0 - M0^(£, V ^ y2 ^) ^ °- 

The Liu equations require, that the entropy must not depend on the gradients 
of the basic state. However, the dissipation inequality still can be solved, if one 
considers two additional physical requirements, prerequisites of relocalizability 

(1) £ is a dynamic variable in a thermodynamic sense therefore £ is zero in 
equilibrium, 

(2) there is no entropy flow connected to the dynamic variable if its value is 
zero. 

It was argued that these requirements are week from a physical point of view. 
With these assumptions one can specify jo in the entropy current with the Nyfri- 
form as j„(£, V£) = A(£, V£)£, or equivalently j„(f , V£) = A(£, V£)d 5 s [TO], accord- 
ing to the mean value theorem. Here the current multipliers A or A are constitutive 
functions to be determined. With the second form of jo the dissipation inequality 
(|13fl can be solved and gives the thermodynamic Ginzburg-Landau equation (JJJ. 
Introducing an additional, new dynamic variable one can recover the additional 
Gurtin-term in the thermodynamic Ginzburg-Landau equation 0|. Let us observe, 
that the key assumption determining the form of the additional entropy current 
was that the entropy current should not affect the equilibrium solutions. Because 
the entropy is a local function in this case, the special internal variables, the current 
multipliers in a sense relocalize the nonlocalities. 

We can apply similar reasoning in the previous non relocalizable case, too. How- 
ever, first we should extend the previous constitutive state space considering a 
derivative one order higher then before. Therefore, let be our constitutive space 
spanned by (£, V£, V 2 £, V 3 £). After a short calculation we can recover the validity 
of in these new variables. The only difference is that the final constitutive 
quantities e.g. jo and L will depend on the larger constitutive state. Now we 
assume that jo has the following form 

jo& V£, V 2 = B(£, V£, V 2 OHt V£, V 2 £), 

were B is a current multiplier. The above form is a direct application of the 
requirement that the entropy current should not change the equilibrium solutions 
(with some minor additional restrictions on the possible constitutive dependencies). 
In this case the entropy production (jl 1|> is 

(14) B • VT+ (V • B + V • d 2 s - dis) ■ T> 0. 

With two undetermined constitutive functions (B, T) the inequality has a general 
solution, the currents and forces are determined by the constitutive dependencies. 
In case of isotropic materials 

B = LrfT 

T = L 2 (V-B + V-d 2 s-dis). 

Here L\ and L 2 are nonnegative scalar constitutive functions. B can be elimi- 
nated from the above equations and finally we get 

(15) dt£ = T = L 2 (V • 8 2 s - d lS ) + L 2 V • (LiVftf)- 

This is the Ginzburg-Landau equation with a characteristic additional Gurtin- 
term. 
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4. Conclusions and discussion 



The requirement of a nonnegative entropy production is a relatively strong and 
not a complete form of the Second Law. Strong form because it is a local require- 
ment and other weaker formulations require only the validity of integral inequalities. 
Not a complete one because an increasing entropy is only a part of the physical con- 
tent of the Second Law. The stability of materials in isolated systems incorporates 
some other conditions (e.g. concave entropy function), too [201 - Colcman-Mizcl 
methodology is a kind of basic philosophical requirement of a thermodynamic the- 
ory: the acceptable theories are those, where the entropy inequality is the conse- 
quence of pure material properties and independent of other elements of the theory 
(e.g. initial conditions) ensuring a kind of universality and some stability properties 
to any thermodynamic theories. 

It is interesting to know, that the doubled variational-thermodynamic structure 
of the Ginzburg-Landau equation can be generalized considerably. That is the 
idea behind the General Equation for the Nonequilibrium Reversible-Irreversible 
Coupling (GENERIC) EI], where the variational part and a formalism 

from mechanics plays the leading role (different brackets, geometrical point of view, 
etc.), but both parts are represented. In this paper we unified the variational and 
the thermodynamic parts of the derivation of the Ginzburg-Landau equation on 
a pure thermodynamic ground, where we did not refer to any kind of variational 
principle. However, the derived static part turned out to have a complete Euler- 
Lagrange form. The dynamic part contains a first order time derivative therefore 
one cannot hope to derive it from a variational principle of Hamiltonian type |24| . 
In our approach we get the "reversible", "variational" part as a specific case of 
the thermodynamic, irreversible thinking, but one cannot hope the contrary, the 
irreversible part cannot be derived from a variational, reversible thinking. 

Weakly nonlocal, pattern forming equations emerge in different fields of physics 
independently of thermodynamic argumentation. Understanding their compatibil- 
ity with the Second Law can be considered as one of the most important challenges 
of contemporary non-equilibrium thermodynamics. In this paper it was shown that 
one of the most important pattern forming equations, the Ginzburg-Landau equa- 
tion, is a straightforward consequence of the entropy inequality alone in a nonlocally 
extended constitutive space. On the other hand the outline of the mathematical 
background and all the key ingredients of an efficient formalism to exploit the 
Second Law inequality in weakly nonlocal continuum theories (gradient theories, 
theories with coarse grained thermodynamic potentials, phase-field models, etc..) 
is given. 

5. Appendix: Liu theorem as a variant of Farkas' lemma and some of 

its consequences 

In 1972 Liu introduced a method of the exploitation of the entropy principle 
|25| . Liu's procedure became a basic tool to find the restrictions posed by the 
entropy inequality. The method is based on a linear algebraic theorem, called 
Liu's theorem in the thermodynamic literature [HI 1 1 -ij and on an interpretation 
of the entropy inequality, one of the fundamental ingredients of the Second Law. 
Recently Hauser and Kirchner recognized that Liu's theorem is a consequence of 
the fundamental theorem of linear inequalities, a famous statement of optimization 
theory and linear programming, the so called Farkas' lemma |2fi| . That theorem 
was proved first by Farkas in 1894 P2| and independently by Minkowski in 1896 
|28j . In this appendix we formulate and generalize Liu's theorem in a way that is 
best adapted for our purposes and shows the whole train of thought from Farkas' 
lemma to Liu's theorem giving a simple proof to every statement in question. 
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Farkas' lemma can be formulated in several different forms, that are more or less 
equivalent |29l \'M)\ . Here we start from a simple variant. 

Lemma 5.1. (Farkas) Let a^ ^ be independent vectors in a finite dimensional 
vector space V, i — l...n, and S — {p G V*|p • a^ > 0, i = l...n}. The following 
statements are equivalent for a b G V: 

(i) P • b > 0, for all p G S. 

(ii) There are nonnegative real numbers X%, A„ such that b = Y^i=i -^ a i- 

Proof: S is not empty. In fact, for all k,i G {1, ...,n} there is a p& G V* such 
that p*; • afe = 1 and pj, ■ a; = if i ^ fc. Evidently pj, G S for all k. 

(ii) =*> (i) p • ^i a i = S™=1 A iP ■ a; > if p G S". 

(i) (ii) Let 5 = {y G V*|y • a, = 0,i = l...n}. Clearly ^ S C 5. 

If y G So then — y is also in So, therefore y • b > and y • b > together. 
Therefore for all y G So it is true that y • b = 0. 

As a consequence b is in the set generated by {a^}, that is there are real numbers 
Ai, A„ such that b = X)"=i ^i a i- These numbers are nonnegative, because with 
the previously defined pfe G S, < pk ■ b = p/,, ■ Yli=i -^ a « = \Pk ■ a ; = is valid 
for all fc. ■ 

Remark 5.1. In the following the elements ofY* are called independent variables 
and V* itself is called the space of independent variables. The inequality in the 
first statement of the lemma is called aim inequality and the nonnegative numbers 
in the second statement are called Lagrange- Farkas multipliers. The inequalities 
determining S are the constraints. 

In the calculations an excellent reminder is to use Lagrange- Farkas multipliers 
similarly to Lagrange multipliers in case of conditional extremum problems: 

n / n \ 

pb-^A i pa i =p b-^A i a i >0, Vp G V* 

i=l \ i=l / 

From this form we can read out the second statement of the lemma. 

Remark 5.2. The original statement does not require the independency of the 
vectors in the constraint. We need some extra conditions and that generalization 
destroys the simplicity of the proof. However, we do not need this generalization in 
thermodynamics. 

The geometric interpretation of the theorem is important and graphic: either the 
vector b belongs to the cone generated finitely by the vectors a^ (Cone(a.\, ...,a n ) = 
{Aiai + ... + A n a n | (Ai, A„) G R +n ), or there exists a hyperplane separating b 
from the cone. 

5.1. Affine Farkas' lemma. This generalization of the previous lemma was first 
published simultaneously by A. Haar and J. Farkas in the same number of the 
same journal, with different proofs \M\ YA2\ . Later it was reproved independently 
by others several times (e.g. |33l !3(Jp . Here we give a simple version again. 

Theorem 5.2. (Affine Farkas) Let a^ ^ be independent vectors in a finite di- 
mensional vector space V and ai real numbers, i — l...n and Sa — {p G V*|p • a^ > 
cti,i = l...n\. The following statements are equivalent for a b G V and a real 
number j3: 

(i) p • b > [3, for all p G Sa- 

(ii) There are nonnegative real numbers Ai, A„ such that h — Y^7=i ^ iSLi an< ^ 

Proof: Sa is not empty. In fact, aiPk G Sa for ah k (pk ■ a*. = 1 and p^ • a,; = 
if i ^ k as previously). 
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(a) (i) p • b = p • Yh=i Xi&i = 127=1 ' &i - 127=1 ^ 

(i) (m) First we will show indirectly that the first condition of lemma l5~TI is a 
consequence of the first condition here, that is if (i) is true then p • b > 0, for all 
peS. 

Thus let us assume the contrary, hence there is p' G S, for which p' • b < 0. 
Take an arbitrary p e Si, then p + kp' G Sa for all real numbers k. But now 
(p + fcp') • b = p • b + kp' ■ b < (3, if k > ■ That is a contradiction. 

Therefore, according to Farkas' lemma (Lemma |5.1(l there exist Lagrange-Farkas 
multipliers A = (Ai, A n ) G M™ + such that b = Y17=i Aj a i- Hence (3 < inf p& s A {p- 

127=i Ai*} = in fp£s A {127=i a»p • a 'i = 127=1 ^ a i- m 

Remark 5.3. The multiplier form is a good reminder again 

n n n 

(p-b-(3)-J2 MP • »i - on) = P ■ (b - X * ■ a ~ P + J2 Xi0li ^ °' Vp e V * 

i— 1 i— 1 i— 1 

Remark 5.4. The geometric interpretation is similar to the previous one, but 
everything is affine. 

5.2. Liu's theorem. Here the constraints are equalities instead of inequalities, 
therefore the multipliers are not necessarily positive. 

Theorem 5.3. (Liu) Let a^ ^ be independent vectors in a finite dimensional 
vector space V and ai real numbers, i = l...n and Sl = {p £ V*|p ■ * = on,i = 
l...n}. The following statements are equivalent for a b G V and a real number (3: 

(i) p • b > [3, for all p G Sl, 

(ii) There are real numbers Ai, A„ such that 

n 

(16) b = A,-* 

i=l 

and 

n 

(17) f3<J2 X ^- 

Proof: A straightforward consequence of the previous affinc form of Farkas' 
lemma because Sl can be given in a form Sa with the vectors a^ and — a;, i — 
1, ...,n: S L = {p G V*|p • a t > a it and p • (-a*) > -a,, i = l...n}. 

Therefore there are nonnegative real numbers , A+ and Aj~, A~ such, that 

b = £™ =1 (a^* - a,"*) = Er=i(A+ - K)°* = 127=1 A * a * and f 3 ^ 127=i( x t a i - 

Remark 5.5. The multiplier form is a help in the applications again 

n n n 

< (p-b-^)-^A i (p-a i -a i ) = p-(b-^A i -a i )- ( 9 + ^A i a i! Vp G V* . 

i=l i=l i=l 

Remark 5.6. Ln the theorem with Lagrange multipliers for a local conditional ex- 
tremum of a differentiable function we apply exactly the above theorem of linear 
algebra after a linearization of the corresponding functions at the extremum point. 

Considering the requirements of the applications we generalize Liu's theorem 
to take into account vectorial constraints. First of all let us remember some well 
known identifications of linear algebra: Lin(W , V) = Bilin(W x V*, R) = V <X> U, 
where Bilin denotes the bilinear mappings of the corresponding spaces (see e.g. 

EH)- 
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Theorem 5.4. (vector Liu) Let A ^ in a tensor product V ® U of finite dimen- 
sional vector spaces V and U. Let a € U and SVl = {p £ V* |p • A = a}. The 
following statements are equivalent for a b G V and a real number (3: 

(i) p • b > (3, for all p G SVl. 

(zi,) There is a A G U* suc/i i/iai 

(18) b = A A, 
and 

(19) [3<\a. 

Proof: Let us observe that we can get back the previous form of the theorem 
by introducing a linear bijection K : U — > K n , a coordinatization in U. Therefore, 
applying it for K • A = (A)j = a^, K • a. — (a)i = on and K' • A = K' • a. — a[ 
we get that b = £™ =1 A i&i = £* =1 A^. Thus A- = K'*" 1 K* ■ A;. Therefore there 
is a A G U, independently of the coordinatization, with the components Ai and AJ 
in the coordinatizations K and K'. ■ 

The previously excluded degenerate case of A = deserves a special attention. 
Now we require the validity of the aim inequality for all p G V* without any 
constraint. The consequences can be formulated as previously and the proof is 
trivial. 

Theorem 5.5. (degenerate Liu) The following statements are equivalent for a b G 
V and a real number (3: 

(i) p-h> f3 for all p G V*. 

(ii) b = and (3 < 0. 

Remark 5.7. The practical application rule is that if A = then the multiplier is 
zero. 

Remark 5.8. In continuum physics and thermodynamics the corresponding form of 
IJiSj) and \1$) are called Liu equation(s) and the dissipation inequality, respectively. 
We apply the same names for the degenerate case, too. There the Lagrange- Farkas 
multipliers are called simply Lagrange multipliers. Our nomenclature honors Farkas 
and emphasizes the difference between the two kind of multipliers. It can be impor- 
tant also to make a clear distinction of a similar but different nomenclature and 
method in variational principle construction in continuum physics . 
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